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At a critical point of a second order phase transition the intrinsic energy surface 
is flat and there is no stable minimum value of the deformation. However, for a 
finite system, we show that there is an effective deformation which can describe 
the dynamics at the critical point. This effective deformation is determined by 
minimizing the energy surface after projection onto the appropriate symmetries. We 
derive analytic expressions for energies and quadrupole rates which provide good 
estimates for these observables at the critical point. 
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Dynamical systems such as nuclei can undergo phase transitions associated with a change 
of shape of their equilibrium configuration. It has been recognized recently [1, 2] that it is 
possible to formulate the concept of a dynamical symmetry (here meaning solvability in terms 
of quantum numbers) for systems at the critical point of such shape phase transitions. The 
importance of these critical-point symmetries lies in the fact that they provide a classification 
of states and analytic expressions for observables in regions where the structure changes 
most rapidly. For nuclei, two critical-point symmetries, called E(5) [1] (not the Euclidean 
group in five dimensions), and X(5) [2], were considered in the geometric framework of the 
collective model. This model involves a Bohr Hamiltonian which describes the dynamics of 
a macroscopic quadrupole shape via a differential equation in the intrinsic quadrupole shape 
variables (3 and 7. The E(5) [X(5)] critical-point symmetry (CPS) is applicable to a second- 
[first-] order shape phase transition between spherical and deformed 7-unstable [axially- 
symmetric] nuclei. In the present work we focus on the E(5) CPS for which an empirical 
example has been found in 134 Ba [3, 4] and possibly in 104 Ru [5], 102 Pd [6] and 108 Pd [7]. To 
apply the E(5) CPS to real nuclei one has to take into account the finite number of nucleons. 
This can be conveniently done in the algebraic framework of the interacting boson model 
(IBM) [8]. This model describes low- lying quadrupole collective states in nuclei in terms of a 
system of N monopole (s) and quadrupole (d) bosons representing valence nucleon pairs. The 
IBM Hamiltonian relevant to the critical point of the phase transition between spherical and 
7-unstable deformed nuclei preserves the 0(5) symmetry [9]. Its energy surface, obtained 
by the method of coherent states [9, 10], is 7-independent and exhibits a flat-bottomed 
behavior in (3 which resembles the infinite square-well potential used to derive the E(5) CPS 
in the geometric approach [1]. Calculations with finite N values (iV=5 for 134 Ba) have found 
that this critical IBM Hamiltonian can replicate numerically the E(5) CPS and its analytic 
predictions. [3, 5, 6, 7]. In the present work we examine the properties and conditions 



that enable features of E(5) CPS to occur in a finite system. For that purpose we propose 
wave functions of a particular analytic form, which can simulate accurately the exact IBM 
eigenstates at the critical point. These wave functions with fixed N and good 0(5) symmetry 
are used to derive accurate estimates for energies and quadrupole rates at the critical point 
without invoking large- N approximations. The proposed wave functions can be obtained by 
projection from intrinsic states with an effective /^-deformation. 

In the geometric approach the E(5) eigenfunctions [1] are proportional to Bessel functions 
of order r + § and the corresponding eigenvalues are proportional to (x^ T ) 2 . Here r is the 
0(5) quantum number, and x^ T is the £-th root of these Bessel functions. A portion of an 
E(5)-like spectrum is shown in Fig. 1. It consists of states, L^ T , arranged in major families 
labeled by £ = 1, 2, . . . and 0(5) r-multiplets (r = 0, 1, . . .) within each family. The angular 
momenta L for each r-multiplet are obtained by the usual 0(5) D 0(3) reduction [11]. The 
E(5) CPS leads to analytic parameter-free predictions for energy ratios and B(E2) ratios 
which persist when carried over to a finite-depth potential [12]. As seen in Table I, the 
E(5) predicted values are in-between the values expected of a spherical vibrator [U (5)] and 
a deformed 7-unstable rotor [0(6)]. 

In the algebraic approach, the U(5)-0(6) transition region is modeled by the Hamiltonian 

H = eh d + ^A [dt.<2t_( s t)2] [h.c] (1) 

with e and A positive parameters. Here h d is the ci-boson number operator, H.c. stands 
for Hermitian conjugate and the dot implies a scalar product. In the £7(5) limit (A = 0), 
the spectrum of H is harmonic, en d , with n d = 0, 1, 2 ... N. The eigenstates are classified 
according to the chain U(<6) D U(5) D 0(5) D 0(3) with quantum numbers | N,n d ,r, L) 
(for r > 6 an additional multiplicity index is required for complete classification). These 
states can be organized into sets characterized by n d = r + 2k. States in the lowest-energy 
set (k = 0) satisfy P \N,n d = r,r,L) = with P f = S ■ d) . Other sets (k > 0) are 



generated by \N,n d ,r,L) oc (P f ) fe | N - 2k, n d = t,t,L). In the 0(6) limit (e = 0), the 
spectrum is \A(N - a) ( iV + a + 4) with a = N, N - 2, TV - 4, ... or 1. The eigenstates 
are classified according to the chain £7(6) D 0(6) D 0(5) D 0(3) with quantum numbers 
| iV, a, r, L). The ground band has a = N and its members satisfy Pi \ N, a = N, r, L) = 
with P\ = [S ■ S - (st) 2 ]. The remaining bands with a = N — 2k are generated by 
| N,a,r,L) oc (Pi) k | iV — 2k, a, r, L). These results suggest that in-between the £7(5) and 
0(6) limits, we consider a ground band (£ = 1) for the Hamiltonian (1) determined by the 
condition 



P y \£ = l;y,N,T,L) = , 



pt = \s.s-y(jy 



(2) 



In the £7(5) basis these states are 



( = 1; N, r, L) = ]T - [l + (-1)"—] £„ d , T | iV, Wd , r, L> , 



(3) 



and the summation covers the range r < n d < N '. The coefficients £„ d)T have the explicit 
form 

1/2 



(7V-r)!(2r + 3)!! 



(N - n d )\ (n d - t)\\ (n d 



3)! 



(n d -r)/2 



(2r + 3)!! 







a \ 2p (2p+l)!! 



Members of the first excited band (£ = 2) have approximate wave functions of the form 



(4) 



t = 2;y,N,T,L)=MpPl\t = l;y,N-2,T,L) 



Mp=\2(2N + y 2 + l)+A{y 2 -l)S { l 



(N-2) 

T 



-1/2 



(5) 



where S[ J is defined in Eq. (7) below. 



The states of Eqs. (3) and (5) have fixed N, L and good 0(5) symmetry r. Henceforth, for 
reasons to be explained below, they will be referred to as r-projected states. Diagonal matrix 
elements of the Hamiltonian (1) in these states, denoted by E^ T = (£;y, N,r\H\£;y, N,t), 
can be evaluated in closed form 



e[N-SW] + \A(l-yfsV 

e {N - 2Nj [2y* + (2N + 7y> - 1) ^~ 2) + 2(y 2 - 1) ^~ 2) ] } 
+\a2N} { 2(y - If (y 2 - 1) S^- 2) + (y~ I) 2 (2N + y* - 8y + 5) Sg. 

+16(|/ -1)(N + y) S[ N - 2) + 2 [ (2N + y){2N + y + 2) + 1 ]} . (6) 



2 n(N) 



-2) 



The E^ T are independent of L since H is an 0(5) scalar. Their expressions involve the 
quantities Sj^ — (£ — 1; y, AT, r|(st) fc s fc |£ = 1; y, N, r) which are given by 



g^(y) 
G&2n(y) 



(7) 



with Sg° = 5{y5{JJ _1) and sg° = S{ N T ] S{ N T ~ 2) . Non-diagonal matrix elements of H 
between r-projected states in different £-bands, H\,i-, T — (£ — 2;y, N,t\H\£ = l;y,N,r), 
can be evaluated as well 



#1,2;. = 2A/; 



q (N) 



1/2 



1 



{ey + ^A(i,- 1) (2 A' - - y f I ) + 2(,y - 1) 5', 



(7V-2) 



(8) 



By techniques similar to that employed in the 0(6) limit of the IBM [11], explicit expressions 
of quadrupole rates can be derived for transitions between the r-projected states. For the 
relevant IBM quadrupole operator, T(E2) = efts + s^d, these transitions are subject to the 
0(5) selection rule Ar = ±1, and, as explained in the caption of Fig. 1, it is sufficient to 
focus on the B(E2) values of the type 



B{E2- £ = 1; r + 1, L = 2r + 2 -> £ = 1, r, L = 2r) = 



(r + 1) 



tl (q(N)\ 2 G N-r+ 

-T+l) WW 



(r+1) 
-r+1 



(2r + 5)(iV 



(*) 



y + (N-r)^M 



n 2 



B(E2; f = 2, r, L = 2r -> £ = 1, r + 1, L = 2r + 2) = 

(r+l)(4r + 5) giggly) , Q . 

(4r + l)(2r + 5) 4 ^^ ( ^ " 1} (iV " T) ' (9) 

The states in Eq. (3) and Eq. (5) can be obtained by 0(5) projection from the IBM 
intrinsic states for the ground band 



c-N) = (N\)-^ 2 (bl) N \0] 



cos 7 dl + (3 sin 7 — j= (d\ + dL 2 ) + 



(10) 



bt = (i + pr i/2 

and respectively for the (3 band 

\(3;N) = NpPl\c-N-2) , (11) 

provided y = f3 2 . The expressions in Eqs. (3)-(9) depend on the so far unspecified parameter 
y. Normally, the equilibrium value of (3, and hence y, is chosen as the global minimum of 
the intrinsic energy surface determined from the expectation value of H in the intrinsic state 
(10). This is a standard procedure for a Hamiltonian describing nuclei with rigid shapes, 
for which the global minimum is deep and well-localized. However, near the critical point 
of the phase transition an alternative procedure is required. 

The IBM Hamiltonian, H cri , at the critical point of the U(5)-0(6) phase transition cor- 
responds [9] to a special choice of parameters in the Hamiltonian of Eq. (1) 

H cri : e — (N — I) A . (12) 

The intrinsic energy surface of H cri has the form 

E((3) = E + A N(N - 1) + f3 2 )~ 2 (13) 

with E = ^AN(N — 1) a constant. The energy surface E{(3), shown in Fig. (2a), has a 
flat behavior (~ (3 A ) for small /3, an inflection point at f3 — 1 and approaches a constant for 
large (3. The global minimum at (3 = is not well-localized and E(f3) exhibits considerable 
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instability in f3, resembling a square- well potential for < (3 < 1. Under such circumstances 
fluctuations in (3 are large and play a significant role in the dynamics. Some of their effect can 
be taken into account by introducing into the intrinsic states of Eqs. (10) and (11) an effective 
/9-deformation. The effective deformation is expected to be in the range < y = (3 2 < 1, 
in-between the respective U(5) and 0(6) value of (3. This will enable a reproduction of E(5) 
characteristic signatures which are in-between these limits (see Table I). In contrast to E(0), 
we see from Fig. (2b) that the 0(5) projected energy surface E^ = i jT=0 (y) of H cri (Eqs. (6) 
and (12) with N = 5), does have a stable minimum at a certain value of y, which we interpret 
as an effective /5-deformation. This procedure, based on variation after projection, is in the 
spirit of [13] in which it is shown that in finite boson systems, a 7-unstable 0(6) state can 
be generated from a rigid triaxial intrinsic state with an effective 7-deformation of 30°. In 
the present case the 7-instability is treated exactly by means of 0(5) symmetry, while the 
/^-instability is treated by means of an effective deformation. The appropriate value of y can 
be used to evaluate the band- mixing, rj T (y) = Jf 1 ^^ ■ ^ small value of rj T will ensure that 
the r-projected states of Eqs. (3), (5) form a good representation of the actual eigenstates 
of H^, an d turn the expressions of Eqs. (6), (9) into meaningful estimates for energies and 
quadrupole transition rates at the critical point. 

To test the suggested procedure we compare in Table II the U(5) decomposition of exact 
eigenstates obtained from numerical diagonalization of H cri for N = 5 with that calculated 
from the r-projected states with y = 0.314 [the global minimum of E 1:0 (y)]. As can be seen, 
the latter provide a good approximation to the exact eigenstates (the corresponding band- 
mixing is 7] T = 0.12, 3.53, 4.14, 3.05% for r = 0, 1, 2, 3). This agreement in the structure of 
wave functions is translated also into an agreement in energies and B(E2) values as shown 
in Table I. The results of Table I and II clearly demonstrate the ability of the suggested 
procedure to provide analytic and accurate estimates to the exact finite-N calculations of 
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the critical IBM Hamiltonian, which in-turn agree with the experimental data in 134 Ba and 
captures the essential features of the E(5) critical-point symmetry. 

To summarize, in this work we have considered properties of a critical-point symmetry 
in a finite system. We have focused on the E(5) critical-point symmetry corresponding to 
flat-bottomed potentials as encountered in a second-order phase transition between spher- 
ical and deformed 7-unstable nuclei. We have shown that intrinsic states with an effective 
/3-deformation, reproduce the dynamics of the underlying non-rigid shapes. The effective 
deformation can be determined from the global minimum of the energy surface after projec- 
tion onto the appropriate symmetry. In the present case, states of fixed N and good 0(5) 
symmetry projected from these intrinsic states provide good analytic estimates to the exact 
eigenstates, energies and quadrupole transition rates at the critical point. 

This work was supported in part by the Israel Science Foundation (A.L.) and in part by 
the U.S. Department of Energy under contract W-7405-ENG-36 (J.N.G). 
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TABLE I: Excitation energies (normalized to the energy of the first excited state) and B(E2) 
values (in units of B(E2; 2+^ — > 0+ ) = 1) for the E(5) critical-point symmetry [1], for several 
N=5 calculations and for the experimental data of 134 Ba [14]. The finite-N calculations involve the 
exact diagonalization of the critical IBM Hamiltonian (H cr i) [Eq. (12)], r-projected states for H cr i 
[Eqs. (6), (9) with y = 0.314], the 17(5) limit [en d ] and the 0(6) limit [(A/4)(N - a)(N + a + 4) + 
Br(r + 3)]. 
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TABLE II: U(5) decomposition (in %) of the L^ T states for N = 5. The calculated values are 
obtained from the r-projected states, Eqs. (4), (5) with y = 0.314. The exact values are obtained 
from numerical diagonalization of the critical IBM Hamiltonian H cr i, Eq. (12). 
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FIG. 1: An E(5)-like spectrum of states labeled by L^ jT . Shown are the transitions whose B(E2) 
values are given in Eq. (9). The E2 rates for other At = 1 transitions (not shown) are governed 
by 0(5) symmetry. Specifically B(E2; h\ 2 -► 2^) for L = 4,2 are in the ratio 1 : 1 respectively 
B(E2; L+3 -► 4+ 2 ) for L = 6,4,3 and B(E2; L+ 3 -► 2+ 2 ) for L = 4, 3, are in the ratios 1 : 
10/21 : 2/7 : 11/21 : 5/7 : 1, respectively. 

FIG. 2: Energy surfaces of the critical IBM Hamiltonian H cr i (12) with N = 5 and A = 1. (a) In- 
trinsic energy surface E(j3), Eq. (13) [solid line], and its approximation by a square- well potential 
[dashed line], (b) 0(5) projected energy surface ^ =ljT=0 (y), Eq. (6). The global minimum is at 
y = 0.314. 



